Spinor Casimir Densities for a Spherical Shell in the 
Global Monopole Spacetime 



A. A. Saharian^'^* and E. R. Bezerra de Mello^^ 
^Department of Physics, Yerevan State University, 
375049 Yerevan, Armenia 
'^the Abdus Salam International Centre for Therietical Physics, 

34014 Trieste, Italy, 
^ Departamento de Fisica-CCEN, Universidade Federal da Paraiha 
58.059-970, J. Pessoa, PB C. Postal 5.008, Brazil 

February 1, 2008 



Abstract 

We investigate the vacuum expectation values of the energy-momentum tensor 
and the fermionic condensate associated with a massive spinor field obeying the 
MIT bag boundary condition on a spherical shell in the global monopole spacetime. 
In order to do that it was used the generalized Abel-Plana summation formula. As 
we shall see, this procedure allows to extract from the vacuum expectation values 
the contribution coming from to the unbounded spacetime and explicitly to present 
the boundary induced parts. As to the boundary induced contribution, two distinct 
situations are examined: the vacuum average effect inside and outside the spherical 
shell. The asymptotic behavior of the vacuum densities is investigated near the 
sphere center and surface, and at large distances from the sphere. In the limit of 
strong gravitational field corresponding to small values of the parameter describing 
the solid angle deficit in global monopole geometry, the sphere-induced expectation 
values are exponentially suppressed. As a special case we discuss the fermionic 
vacuum densities for the spherical shell on background of the Minkowski spacetime. 
Previous approaches to this problem within the framework of the QCD bag models 
have been global and our calculation is a local extension of these contributions. 
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1 Introduction 



Different types of topological defects [1] may have been formed during the phase transi- 
tions in the early universe. Depending on the topology of the vacuum manifold Ai these 
are domain walls, strings, monopoles and textures corresponding to the homotopy groups 
TTo{M), TTi{M), TT2{M) and 7rs{M), respectively. Physically these topological defects ap- 
pear as a consequence of spontaneous breakdown of local or global gauge symmetries of the 
system composed by self-coupling scalar Higgs or Goldstonc fields, respectively. Global 
monopoles are spherically symmetric topological defects created due to phase transition 
when a global symmetry is spontaneously broken and they have important role in the 
cosmology and astrophysics. 

The simplest theoretical model which provides global monopoles has been proposed a 
few years ago by Barriola and Vilenkin [2]. This model is composed by a self-coupling iso- 
scalar Goldstonc field triplet 0", whose original global 0(3) symmetry is spontaneously 
broken to U{1). The matter field plays the role of an order parameter which outside the 
monopole's core acquires a non- vanishing value. The main part of the monopole's energy 
is concentrated into its small core. Coupling this system with the Einstein equations, a 
spherically symmetric metric tensor is found. Neglecting the small size of the monopole's 
core, this tensor can be approximately given by the line element 

ds' = dt" - dr" - aV^ ^ sin^ ^#2^ , (1) 

where the parameter a^, smaller than unity, depends on the symmetry breaking energy 
scale and codifies the presence of the global monopole^. This spacetime corresponds to an 
ideahzed point-like global monopole. It is not flat: the scalar curvature R — 2(1 — Q;^)/r^, 
and the solid angle of a sphere of unit radius is Q = Ana"^, so smaller than the ordinary 
one. The energy-momentum tensor associated with this object has a diagonal form and 
its non-vanishing components read Tq = = (a^ — l)/r^. 

The quantum effects due to the point-like global monopole spacetime on the matter 
fields have been considered in Refs. [3] and [4] to massless scalar and fermionic fields, 
respectively. In order to do that, the scalar and spinor Green functions in this background 
have been obtained. More recently the effect of the temperature on these polarization 
effects has been analysed in [5] for scalar and fermionic fields. The calculation of quantum 
effects on massless scalar field in a higher dimensional global monopole spacetime has also 
been developed in [6]. 

Although the deficit solid angle and also curvature associated with this manifold pro- 
duce non- vanishing vacuum polarization effects on matter fields, the influence of boundary 
conditions obeyed by the matter fields on the vacuum polarization effects have been inves- 
tigated. The Gasimir energy associated with massive scalar field inside a spherical region 
in the global monopole background have been analyzed in Refs. [7, 8] using the zeta func- 
tion regularization method. More recently the Gasimir densities induced by a single and 
two concentric spherical shells have been calculated [9, 10] to higher dimensional global 
monopole spacetime by making use of the generalized Abel-Plana summation formula 
[11, 12]. This procedure allows to develop the summation over all discrete mode. Here 
we shall calculate the Gasimir densities for fermionic fields obeying MIT bag boundary 

^In fact the parameter = 1 — SirGrj'^, being 77 the energy scale where the global symmetry is 
spontaneously broken. 
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condition on the spherical shell in the point-like global monopole spacetime. Specifically 
we shall calculate the renormalized vacuum expectation values of the energy-momentum 
tensor and the fermionic condensate in the regions inside and outside the spherical shell. 
As we shall see using the generalized Abcl-Plana summation formula, all the compo- 
nents of the vacuum average of the energy-momentum tensor can be separated in two 
contributions: boundary dependent and independent ones. The boundary independent 
contribution is similar to previous result obtained in [4] using different approach. It is 
divergent and consequently in order to obtain a finite and well defined expression we must 
apply some regularization procedure. The boundary dependent contribution is finite at 
any strictly interior or exterior point and docs not contain anomalies. Consequently, it 
does not require any regularization procedure. Because the analysis of boundary inde- 
pendent term has been performed before, in this present analysis we shall concentrate on 
the boundary dependent part. Taking a — 1, from our results in this paper as a spe- 
cial case wc obtain the fermionic Casimir densities for a spherical shell on background of 
the Minkowski spacetime. Motivated by the MIT bag model in QCD, the corresponding 
Casimir effect was considered in a number of papers [13, 14, 15, 16, 17, 18, 19, 20] (for 
reviews and additional references see [21, 22, 23, 24]). To our knowledge, the most of 
the previous studies were focused on global quantities, such as the total vacuum energy 
and stress on the surface. The density of the fermionic vacuum condensate for a massless 
spinor field inside the bag was investigated in Ref. [15] (see also [23]). In the considera- 
tions of the Casimir effect it is of physical interest to calculate not only the total energy 
but also the local characteristics of the vacuum, such as the energy-momentum tensor 
and vacuum condensates. In addition to describing the physical structure of the quantum 
field at a given point, the energy- momentum tensor acts as the source of gravity in the 
Einstein equations ^. It therefore plays an important role in modelling a self-consistent 
dynamics involving the gravitational field [25]. For the case of the Minkowski bulk, our 
calculation is a local extension of the previous contributions on the fermionic Casimir 
effect for a spherical shell. 

This paper is organized as follows: In section 2 we obtain the normalized eigenfuntions 
for a massive spinor field on the global monopole spacetime inside a spherical shell of finite 
radius. In section 3, using the generalized Abel-Plana summation formula, we formally 
obtain the vacuum expectation value of the energy-momentum tensor considering that 
the fermionic field obeys the MIT bag condition on the spherical shell. Explicit behavior 
for boundary dependent term is exhibited. The section 4 is devoted to the calculation of 
the vacuum expectation values for the region outside the shell. In section 5 we present 
our concluding remarks and leave for the Appendix some relevant calculations. 

2 Eigenfunctions for a Spinor Field on the Global 
Monopole Spacetime 

The dynamics of a massive spinor field on a curved spacetime is described by the Dirac 
differential equation 

i7'^(9^ + r^)V'-MV' = 0, (2) 

^Thc effects of the back-reaction corrections on the Einstein equation due to the vacuum polarization 
produced by massless scalar field in a global monopole spacetime has been analyzed in [3] . 
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where 'y^ are the Dirac matrices defined in such curved spacetime, and the spin con- 
nection defined as 

= , (3) 

being the standard covariant derivative operator. Notice that, using the usual com- 
mutation relations for the Dirac matrices, {7^, Y} = 2^'^'', we can see that 

= + (9^7. - 9.7m) ■ (4) 

After this brief introduction, let us now specialize in the spacetime associated with 
the point-like global monopole whose line element is described by (1). In order to develop 
such procedure we shall adopt the following representation for the Dirac matrices 

n _ / 1 \ . f a'' 



^-[0 J' 

given in terms of the curved space Pauli 2x2 matrices a'': 

I _ f cos 9 e~^'^ sin 6* \ 2 _ 1 / ~ sin 6* e~^'^ cos 9 
^ ~ \^ e^-^sin^ -cos^ ) ' ^ ~ '^\e^^cos9 sm9 

3 i / -e-^'t' \ 

^ = -T-r-a \ .i4> n " (6) 



arsin^ V ^ ^ , 
These matrices satisfy the relation 

a'a'^-V' + ^^TmpaP, (7) 

where 7'*^ = —g^^ are the spatial component of metric tensor and 7 is the corresponding 
determinant, e^^"^ is the totally anti-symmetric symbol with £^^^ = 1. Here and below 
the latin indices i,k,--- run over values 1, 2, 3. It can be easily checked that with 
these representations the Dirac matrices satisfy the standard anticommutation relations. 
Substituting these matrices into formula (4), we can see that 

rr, = ^7^ ■ (8) 

ar 

Let us write the four-components spinor field ip in terms of two-components ones as 



Assuming the time dependence in the form e"""^*, from (2) one finds the equations for 
these spinors: 

a'^dk^p + ^^(^ • = ^(^ + M)x , (10a) 

<7'4x+^(^-^)x = i{u;-M)^, (10b) 
ar 
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where a = (cr^, cr^, cr^), and h = f/r. The angular parts of the spinors are the standard 
spinor spherical harmonics ^jim whose explicit form is given in Ref. [26]: 



(11) 



where j specifies the value of the total angular momentum, and m its projection, I 
j ± 1/2, I' = 2j — I. Using the formula 



^ji'm — i {n ■ a)Vtjim , 



it can be seen that 



where we use the notation 



ar 



g\r) + - — -g{r) 



ar 



'jlm : 



K 



-(Z + 1), j ^1 + 1/2 
Z, j ^1-1/2 ■ 



(12) 

(13a) 
(13b) 

(14) 



By taking into account these relations, from (10) we obtain the following set of differ- 
ential equations for the radial functions 



f^r) + ^/(r) + M)g{T) 



ar 



g'^r) + ^^g{r) + (a; - M)/(r) 



, 
. 



ar 



They lead to the second order differential equations for the separate functions: 



r (r) + -f{r) + 
r 

g"{r) + -g\r) + 
r 



, , nin + a) 



K[K — a] 



fir) 

g{r) 



, A; = Vu^ 



0, 



M2 



with the solutions 



y(^) ^ ^^ Zi,/a+i/2\ikr) ^ ^^^^ ^ ^^|./a-i/2|(/cr) 



(15a) 
(15b) 

(16a) 
(16b) 

(17) 



where Zy{x) represents the cylindrical Bessel function of the order v. The constants A 
and B are related by equations (15): 



B 



^kA 



(18) 



for 7 = Z ± - 

As a result for a given j we have two types of eigenfunctions with different parities 
corresponding to j = Z ± 1/2. These functions are specified by the set of quantum 
numbers /3 = {akjlm) and have the form 

Ae-"^' ( Z,^{kr)nji, 



I 



j - ^, u)^±E, Vk^ + M2 , 



(19) 
(20) 
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where j = 1/2, 3/2, . . ., and m = —j, ■ ■ ■ ,j, 

<7 = 0, 1, n„ = {-ir, u^^l±ll^-^. (21) 

a 2 

On the base of formula (19) we define the positive and negative frequency eigenfunctions 
as 



= { 4'!' " ^ °' (22) 

^ ^ i^y, for w<0. ^ ^ 



These functions are orthonormahzed by the condition 

J d'x^^^^^^^P^^:^ = 6,,^60p, ,ri,ri' = ±, (23) 
from which the normahzation constant A can be determined. 

3 Vacuum Expectation Values of the Energy- 
Momentum Tensor Inside a Spherical Shell 

In this section we shall consider the vacuum expectation values of the energy-momentum 
tensor inside a spherical shell concentric with the global monopole. The integration in 
formula (23) goes over the interior region of the sphere and Z^{x) = Ju{x), where Ju{x) 
is the Bessel function of the first kind. We shall assume that on the sphere surface the 
field satisfies bag boundary conditions: 

{l + iYn,^)^^0 , r^a, (24) 

where a is the sphere radius, n;,^ is the outward-pointing normal to the boundary. For 
the sphere nb/^ — (0, 1, 0, 0). In terms of the spinors (p and x this condition is written in 
the form 

(/? + i (a ■ nfe) X = , r = a. (25) 

The imposition of this boundary condition on the eigenfunctions (19) leads to the following 
equations for the eigenvalues 

Jya+nA^a) = n„ ^^^^ Ju^ka) . (26) 
This boundary condition can be written in the form 

J..(A)a) = 0, (27) 
where now and below for a given function F{z) we shall use the notation 

F{z) = zF'{z) + (/X + s^^z^ + n''- (-I)V) F{z) , a = 0, 1, (28) 

with — sgn(a;) and = Ma. Let us denote by \va,s — ka, s — 1,2, ... , the roots to 
equation (27) in the right half plane, arranged in ascending order. By taking into account 
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Eq. (26) and using the standard integral for the Bessel functions, from condition (23) for 
the normalization coefficient one finds 



with u = ±-^aJJV«^ + 

Now we expand the field operator in terms of the complete set of single-particle states 



i^u — 



2a{uj + M) 



, z^K^,,, (29) 



(-) 

/3 



(30) 



where is the annihilation operator for particles, and 6^ is the creation operator for 
antiparticles. In order to find the vacuum expectation value for the operator of the 
energy-momentum tensor we substitute the expansion (30) and the analog expansion for 

the operator %[) into the corresponding expression for the spinor fields. 



T^,u 1^,^} = ^ ^7(/.V^)t/' - (V(^'0)7^)'0 



(31) 



By making use the standard anticommutation relations for the annihilation and creation 
operators, for the vacuum expectation values one finds the following mode-sum formula 



(0|T,.|0) = 5^T,.{4-)(a:),4-)(x)} , 



(32) 



where |0) is the amplitude for the corresponding vacuum. Since the spacetime is spher- 
ically symmetric and static, the vacuum energy-momentum tensor is diagonal, moreover 
{Tg) = (T?). So in this case we can write: 



(0 |r;| 0) = diag(£, -p, -px, -p^) , 



(33) 



with the energy density s, radial, p, and azimuthal, p±, pressures. As a consequence of 
the continuity equation V,^ (O IT^I O) = 0, these functions are related by the equation 



r^ + 2(p-p±) = 0. 



(34) 



which means that the radial dependence of the radial pressure necessarily leads to the 
anisotropy in the vacuum stresses. 

Substituting eigenfunctions (19) into Eq. (32), the summation over the quantum 
number m can be done by using standard summation formula for the spherical harmonics. 
For the energy-momentum tensor components one finds 



-1 



(2j + 1) E J2^uAK,s)fit [Xu.,s, J..{K,sr/a)] , g = £, p, , 



(T=0,1 S=l 



(35) 
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where we have introduced the notations 



(37) 



ftHz,My)] = 



JM - ^^^Y^My)J^^nM + J'+nAy) 



ft^[z,JM] = ^^^^^i=^.(2/)^.+n.(2/) . (38) 

Note that in (35) we have used the relation between the normahzation coefficient and the 
function Ti,{z) introduced in Appendix A: 

The vacuum expectation values (35) arc divergent and need some regularization procedure. 
To make them finite we can introduce a cutoff function $,,(2;), z = Xy^^s with the cutoff 
parameter 77, which decreases with increasing z and satisfies the condition — > 1, 77 —> 
0. Now to extract the boundary-free parts we apply to the corresponding sums over s 
the summation formula derived in Appendix A. As a function f[z) in this formula we 
take f{z) = fjj})^ [z, J,,^{zr/a)] ^niz). As a result the components of the vacuum energy- 
momentum tensor can be presented in the form 

(l(r) = qm(r) + qb(r), q^e,p,p± , (40) 

where the first term on the right hand side comes from the integral on the left of summation 
formula (78), and the second term comes from the integral on the right of this formula. 
Making use the asymptotic formulae for the Bessel modified functions, it can be seen that 
for r < a the part qb{r) is finite in the limit 1] ^ and, hence, in this part the cutoff can 
be removed. As it has been pointed out in Appendix A, the function fili [z, J,y^{zr/a)] 
satisfies relation (81) and, hence, the part of the integral on the right of formula (78) over 
the interval (0,yU.) vanishes after removing the cutoff. Introducing the notation 

11 

^ = ^1 = - + :t , (41) 

a Z 

explicitly summing over a and transforming from summation over j to summation over 
I = j + 1/2, one receives 



where we use the notation 



(43) 



fi'^[x,My)] = 4iy) + 4-M--My)Ju-i{y), (44) 

y 

f!)^^^[x,My)] = ^J.(?/)J._i(y). (45) 
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Introducing in the expression for qb{r) the modified Bessel functions, after some transfor- 
mations end exphcitly summing over a, we obtain the formula 



1 



oo 



oo 



x^dx W[h{x),K^{x)] 




,/^^2^W[I,ix),I,{x)] 



Fi^)[x,7.,(xr/a)] , 



(46) 



1=1 



with 




(47) 



(48) 




(49) 



W [f{x),g{x)] = [xf{x) + (/. + u)f{x)] [xg'ix) + (/. + iy)g{x)] + {x' - f^')f{x)g{x). (50) 



As we see, the part gm(^) in the vacuum expectation value for the energy- momentum 
tensor do not depend on the radius of the sphere a, whereas the contribution of the terms 
qb{r) tends to zero as a ^ oo (for large a the subintegrand behaves as e^^^''/'^"^)). It follows 
from here that the quantities (42) are the vacuum expectation values for the components 
of the energy-momentum tensor for the unbounded global monopole space: 



where \0m) is the amplitude for the corresponding vacuum. Note that in expressions (42) 
for the corresponding components we have not explicitly written the cutoff function. Pre- 
cisely speaking in this form all terms related with (42) are divergent. The renormalization 
prescriptions adopted to provide a finite and well defined result to them arc the usual ones 
applied for the curved spacetime without boundary [27, 28, 25]. For the specific system 
analyzed here the point-splitting renormalization procedure has been applied in previous 
publication [4]. The part qb{r) in Eq. (40) is induced by the presence of the spherical shell 
and can be termed as the boundary part. As we have seen, the application of the gener- 
alized Abcl-Plana formula allows us to extract from the vacuum expectation value of the 
energy-momentum tensor the contribution due to the boundary-free monopole spacetime 
and to present the boundary-induced part in terms of exponentially convergent integrals 
(for applications of the generalized Abel-Plana formula to a number of Casimir prob- 
lems with various boundary geometries sec Rcfs. [11, 12, 29, 30, 31, 32, 33]). It can be 
easily checked that the both terms on the right of formula (40), qm{r) and ^'^(r), obey 
the continuity equation (34). In addition, as it is seen from expressions (46)-(49), for a 
massless spinor field the boundary-induced part of the vacuum energy-momentum tensor 
is traceless and the trace anomalies are contained only in the purely global monopole part 
without boundaries. 

Having the components of the energy-momentum tensor we can find the correspond- 
ing fermionic condensate {Q\ipip\Q) making use the formula for the trace of the energy- 
momentum tensor, Tj^ — It is presented in the form of the sum 



{^m\T^\^m) = diag(£^, -Prn, -Pm, -P±m) , 



(51) 




(52) 



9 



where the boundary-free part (first summand on the right) and the sphere- induced part 
(second summand on the right) are determined by formulae similar to Eqs. (42) and (46), 
respectively, with replacements 

Si?n':,Uy)] ^ ^[JM + Jl-M], (53) 

with notation (50). Alternatively one could obtain formulae (52)-(54) applying the sum- 
mation formula (78) to the corresponding mode-sum for the fermionic con- 
densate. 

Note that formulae (40), (42), (46) can be obtained by another equivalent way, ap- 
plying a certain first-order differential operator on the corresponding Green function and 
taking the coincidence limit. To construct the Green function we can use the correspond- 
ing mode expansion formula with the eigenfunctions (19). This function is a 4 x 4 matrix 
and the angular parts of the corresponding elements are products of the components for 
the spinor spherical harmonics, ^^7^(6*, 0)^j"m^(^') 0')' where the upper indices n, n' = 1, 2 
numerate the spinor components. These parts are the same as in the boundary-free case 
and coincide with the corresponding functions for the Minkowski bulk. The radial parts 
for the components of the Green function contain the products of the Bessel functions 
in the forms Ju„+rn„{zr /a)Ji,^+r'n^{zr' /a), t,t' = 0, 1, where z = \ua,s- To evaluate the 
sum over s we can apply the summation formula (78). The condition (75) is satisfied if 
r -\- r' + \t — t'\ < 2a. The term with the integral on the left of formula (78) gives the 
Green function for the boundary-free global monopole spacetime, and the term with the 
integral on the right will give the boundary-induced part. 

In the case a = 1 the quantities (42) present the vacuum expectation values for the 
Minkowski spacetime without boundaries. This can be also seen by direct evaluation. For 
example, in the case of the energy density making use the formula Yl'iloi'^^ ~^ ^) •Ji+i/2iy) ~ 
2y/TT, one finds 

-j^ °0 /"OO 




which is precisely the energy density of the Minkowski vacuum for a spinor field. As for 
the Monkowski background the renormalized vacuum energy-momentum tensor vanishes, 
?m('^)ren = 0, the vacuum energy- momentum tensor is purely boundary-induced and the 
corresponding components are given by formulae (46)-(49) with u — I + 1/2. Note that 
the previous investigations on the spinor Casimir effect for a spherical boundary (see, 
for instance, [13, 14, 16, 17, 18, 19, 20, 21, 22, 23, 24] and references therein) consider 
mainly global quantities, such as total vacuum energy. For the case of a massless spinor 
the density of the fermionic condensate {'4>'il>)b is investigated in [15] (see also [23]). The 
corresponding formula derived in Ref. [15] is obtained from (46) with replacement (54) 
in the limit n = 0. In Fig. 1 we have presented the dependence of the Casimir densities, 
a^qi){r) on the rescaled radial coordinate r/a for a massless spinor field on the Minkowski 
bulk. The vacuum energy density and pressures are negative inside the sphere. 



10 



0.1 




Figure 1: Vacuum energy density, a'^e (curve a), azimuthal pressure a'^p± (curve b), and 
radial pressure a'^p (curve c) for a massless spinor as functions on the ratio r/a inside and 
outside a spherical shell in the Minkowski spacetime (a = 1). 

Now we turn to the consideration of various limiting cases of the expressions for the 
sphere-induced vacuum expectation values. In the limit r — > 0, for the boundary parts 
(46) the summands with a given / behave as r^^/"^^, and the leading contributions come 
from the lowest / = 1 terms. Making use standard formulae for the Bessel modified 
functions for small values of the argument, for the sphere- induced parts near the center, 
r <^ a, one finds 



2 

TT "a ' / r \ I X 



2.-4 roo W[h{x),K,{x)] 



" "^^^ " 2a(2 + a)P(i + i) l2^J J, _ ^2 W [l{x), U{x)\ ^ ^'^^ 

where u = 1/a + 1/2 and r(a;) is the gamma function. Hence, at the sphere center the 
boundary parts vanish for the global monopole spacetime (a < 1) and are finite for the 
Minkowski spacetime (a = 1). Note that in the large mass limit, jj, ^ 1, the integrals 
in Eqs. (56), (57) are exponentially suppressed by the factor e~'^^. In the Minkowski 
background case the vacuum stresses for a massless spinor are isotropic at the sphere 
center and after the numerical evaluation of the integral one finds 

Eb 2 f°° e-^'={x'^ + x-e'=smhx)dx 0.00579 



(2^2 + 1) cosh(2a:) - 2x sinh(2a:) - 1 ' 

(58) 

The boundary induced parts of the vacuum energy-momentum tensor components 
diverge at the sphere surface, r — > a. These divergences are well-known in quantum field 
theory with boundaries and are investigated for various types of boundary geometries 
[34, 35, 36]. In order to find the leading terms of the corresponding asymptotic expansion 
in powers of the distance from the sphere surface, we note that in the limit r ^ a the sum 
over / in (46) diverges and, hence, for small 1 — r/a the main contribution comes from 
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the large values of /. Consequently, rescaling the integration variable x ^ vx and making 
use the uniform asymptotic expansions for the modified Bessel functions for large values 
of the order [37] , to the leading order one finds 



Pb{r) 



II + 1/5 



127r2a(a-r)3' 



1 - - ) P±b{r) 

0, ' 



1/5 -2/i 



247r2a2(a-r)2' 



(59) 
(60) 



Notice that the terms in these expansions diverging as the inverse fourth power of the 
distance have cancelled out. This is a consequence of the conformal invariance of the 
massless fermionic field and is in agreement with general conclusions of Ref. [35]. Near 
the sphere surface the energy density is negative for all values of /i, while the vacuum 
pressures are negative for < 0.1 and are positive for // > 0.1. It is of interest to note 
that the leading terms do not depend on the parameter a and, hence, are the same for the 
Minkowski and global monopole bulks. For the latter case due to the divergences, near the 
sphere surface the total vacuum energy-momentum tensor is dominated by the boundary 
induced parts ^^(r). The dependence of these parts on the rescaled radial coordinate rja 
is depicted in Fig. 2 for the case of a massless fermionic field on the global monopole 
background with the solid angle deficit parameter a = 0.5. 



0.1 



0.05 



-0.05 



0.25 0.5 0.75 1 1.25 1.5 1.75 2 



Figure 2: The same as in Fig. 1 for the boundary-induced parts a'^Qbir) on background 
of the global monopole spacetime with a = 0.5. 



Now let us consider the limit a <^ 1 for a fixed value r < a. This limit corresponds 
to strong gravitational fields. In this case from (41) one has u ^ l/a ^ 1, and after 
introducing in (46) a new integration variable y = x/u, we can replace the modified 
Bessel functions by their uniform asymptotic expansions for large values of the order. 
The integral over y can be estimated by making use the Laplace method. The main 
contribution to the sum over / comes from the summands with / = 1 and the boundary 
parts of the vacuum energy-momentum tensor components behave as exp[— 2 ln(a/r)/a] 
with Pb/p±b ~ tt- Hence, for a <^ 1 the boundary- induced vacuum expectation values are 
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exponentially suppressed and the corresponding vacuum stresses are strongly anisotropic. 
Fig. 3 shows that the nonzero mass can essentially change the behavior of the vacuum 
energy-momentum tensor components. In this figure we have depicted the dependence of 
the boundary induced quantities a^g;,(r) on the parameter Ma for the radial coordinate 
r = 0.5a. The left panel corresponds to the sphere in the Minkowski spacetime [a = 1) 
and for the right panel a = 0.5. 
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Figure 3: Boundary-induced vacuum expectation values a^qb{r), q = e,p,p±, as functions 
on n = Ma for r/a = 0.5. The curves a, b, c correspond to the energy density (e), 
azimuthal pressure and radial pressure (p), respectively. For the left panel a — 1 

(Monkowski spacetime) and for the right panel a — 0.5. 



4 Vacuum Expectation Values Outside a Spherical 
Shell 

Now let us consider the expectation values of the energy-momentum tensor in the region 
outside a spherical shell, r > a. The corresponding eigenfunctions have the form (19), 
where now the function Zy{kr) is the hnear combination of the Bessel functions of the 
first and second kinds. The coefficient in this linear combination is determined from the 
boundary condition (25) and one obtains 

Zy{kr) = gu{ka, kr) = Ji,{kr)Y^{ka) — Y^{kr)Ji,{ka) , (61) 

where Yj^[z) is the Bessel function of the second kind, and the functions with tilda are 
defined as (28). Now the spectrum for the quantum number k is continuous and the 
corresponding 5kk' in Eq. (23) is understood as the Dirac delta function 5{k — k'). To 
find the normalization coefficient A from Eq. (23) it is convenient to take P = P' for all 
discrete quantum numbers. As the normalization integral diverges in the limit k — k', the 
main contribution into the integral over radial coordinate comes from large values of r 
when the Bessel functions can be replaced by their asymptotics for large arguments. The 
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resulting integral is taken elementary and for the normalization coefficient we obtain 

A^A A- - + f62l 

^' " 2aMJl{ka)+Yl{ka)] ' 

Substituting the eigenfunctions (19) into the mode-sum formula (32) and taking into 
account Eqs. (61) and (62), we can see that the vacuum energy-momentum tensor has 
the form (33). The diagonal components are determined by formulae 



87rQ;2.. . 

1/2 



where the expressions for faX [x,giy„{x,xr/a)] are obtained from formulae (36)-(38) by 
replacements 

My) 9v{x,y) , Ju+naiv) Ju+na(y)Yu{ka) - Y^+riAy)Mka) ■ (64) 

To find the parts in the vacuum expectation values of the energy-momentum tensor in- 
duced by the presence of the sphere we subtract the corresponding components for the 
monopole bulk without boundaries, given by Eq. (42). In order to evaluate the corre- 
sponding difference we use the relation 

f!.t [x,^,,(x,xr/a)] _ ( ) J - -- y ii4^fi^J \x,Hi^\xr/a)] , (65) 

where H!f\z), s — 1, 2 are the Hankel functions. This allows to present the vacuum 
energy-momentum tensor components in the form (40) with the boundary induced parts 

1 °° f°° T (t) 

"'<'■> = 16^ E + 1) E E / ''-Wi7fit m-r/a)] . (66) 

/o ^— n 1 o— 1 o U llu„ \X ) 



j=l/2 <T=0,1 5=1,2"" ^■^'^T 



On the complex plane x we can rotate the integration contour on the right of this formula 
by the angle tx/2 for s = 1 and by the angle — 7r/2 for s = 2. The integrals over the seg- 
ments {0,iiJ,) and (0, —in) cancel out and after introducing the Bessel modified functions 
one obtains 

qJr) = ^ T I r .^'^^ WMx),K^{x)] 

^^^^^ ^2^2^3,Z.7^ ,/^^^W[K^{x),Mx)f'^ [x,A.(xr/aJJ. [tl) 

Here the expressions for the functions F^^ [x, K^(y)] are obtained from formulae (47)-(49) 
by replacements Iu{y) — K^{y) and ly-iiy) —K^^i{y). With the same replacements, 
from (52), (53), (54) we can obtain formulae for the fermionic condensate in the region 
outside a sphere. As for the interior region, it can be seen that in the limit of strong gravi- 
tational fields, a <^ 1, the boundary induced vacuum expectation values are exponentially 
suppressed by the factor exp[— (2/q;) ln(r/a)], an the corresponding vacuum stresses are 
strongly anisotropic: Pb/p±b ~ ce. 

In the case a = 1 the renormalized values for the boundary-free parts qm{r) vanish 
and from (67) with v = I + \/2 the components of the vacuum energy- momentum tensor 
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are obtained for the region outside a spherical shell on the Minkowski bulk. Previous 
approaches to this problem have been global and our calculation is a local extension of 
these results. Note that for an infinitely thin spherical shell the total vacuum energy 
for a massless spinor, including interior and exterior parts, is positive, E = 0.0204/a 
[16, 19, 20]. In Fig. 1 we have plotted the dependence of the vacuum energy density 
and stresses on the radial coordinate for a massless spinor field outside a sphere on the 
Minkowski bulk. The same graphics for the boundary-induced expectation values (67) on 
background of the global monopole spacetime with a = 0.5 are depicted in Fig. 2. As 
seen form these figures, the energy density and azimuthal pressure are positive outside a 
sphere, and the radial pressure is negative. The latter has the same sign as for the interior 
region. 

For the case of a massless spinor the asymptotic behavior of boundary part (67) at 
large distances from the sphere can be obtained by introducing a new integration variable 
y = xr/a and expanding the subintegrands in terms of a/r. The leading contribution for 
the summands with a given / has an order (a/r)^''"''^ and the main contribution comes from 
the / = 1 term. Evaluating the standard integrals involving the square of the Mac-Donald 
function, the leading terms for the asymptotic expansions over a/r can be presented in 
the form 

a.(r)- ' r(l + l)r(f + |)/, 

''^'^ ~ 2fvraM4-«^)(2 + a)r3(i + i) Kr) ' ^^'^ 

where 

= 4-^^, L = /p, =1. (69) 

■' 3q; + 2' ■'^ 3q; + 2' "'^^ ^ ^ 

As for the interior components, the quantities (67) diverge at the sphere surface r — a. 

Near the surface the dominant contributions come from modes with large / and by making 
use the uniform asymptotic expansions for the Bessel modified functions, the asymptotic 
expansions can be derived in powers of the distance from the sphere. The leading terms 
of these asymptotic expansions are determined by formulae 

127r2a(r-a)3' 

P,{r) -(^-l)px^(^)-- 24XVr'-a)^ - ^^'^ 

Near the sphere the outside energy density is positive for < 0.04 and is negative for 
II > 0.04. Recall that near the sphere the interior energy density is always negative. As 
we see, the leading terms for the radial pressure are the same for the regions outside and 
inside the sphere. For the azimuthal pressure these terms have opposite signs. In the case 
of the massless spinor field the same is true for the energy density. 



5 Concluding Remarks 

In this paper we have analyzed the fermionic Casimir densities induced by a spherical shell 
in a idealized point-like global monopole spacetime. Specifically the renormalized vacuum 

expectation value of the operator energy-momentum tensor has been considered when the 
matter fields obey the MIT bag boundary condition on the shell. Because the boundary 
condition provides a discrete energy spectrum for the matter fields, the summation over 
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the modes can be developed by using the the generahzed Abel-Plana summation formula. 
Moreover, this procedure allows us to extract from the vacuum average the boundary 
dependent part. This part presents, besides the contribution coming from the parameter 
a which characterizes the presence of the global monopole, contributions coming from 
the boundary itself. Two distinct situations have been considered: the calculation of the 
vacuum average in the regions inside and outside of the spherical shell. As to the inside 
region case, we pointed out that all contributions to the vacuum average go to zero as the 
radius of the spherical shell goes to infinity, as it was expected. 

The boundary-induced expectation values for the components of the energy-momentum 
tensor are given by formulae (46) and (67) for interior and exterior regions, respectively. 
The corresponding formulae for the fermionic condensate densities are obtained from these 
expressions with replacement (54) for the interior region and with additional replacement 
^v{v) — ^ Ki,(y), 7j,_i(y) — > —K^_i{y) for the exterior region. These expressions diverge 
in a non-integrable manner as the boundary is approached. The energy density and 
azimuthal pressure vary, to leading order, as the inverse cube of the distance from the 
sphere, and near the sphere the azimuthal pressure has opposite signs for the interior 
and exterior regions. For a massless spinor the same is true for the energy density. The 
radial pressure varies as the inverse square of the distance and near the sphere has the 
same sign for exterior and interior regions. This behavior is clearly seen from Figs. 1 
and 2 where the radial dependence of the vacuum energy density, azimuthal and radial 
pressures are presented for the Minkowski background and global monopole spacetime 
with a — 0.5. The leading terms of the corresponding asymptotic expansions near the 
sphere do not depend on the solid angle deficit parameter and are the same for these 
two cases. Near the sphere the interior energy density is negative for all values of the 
mass, while the exterior energy density is positive for Ma < 0.04 and is negative for 
Ma > 0.04. in order to illustrate the dependence on the mass, in Fig. 3 we have plotted 
the boundary-induced vacuum densities at r = 0.5a as functions on Ma. Near the sphere 
center the dominant contributions come from modes with I — 1 and the sphere-induced 
vacuum expectation values vanish for the global monopole spacetime and are finite for 
the Minkowski bulk. The asymptotic behavior of the vacuum energy-momentum tensor 
at large distances from the sphere is described by formula (68) and the corresponding 
diagonal components go to zero as (a/r)^/""*"^. In the hmit of strong gravitational field, 
corresponding to small values of the parameter a, describing the solid angle deficit, the 
boundary- induced part of the vacuum energy- momentum tensor is strongly suppressed 
by the factor exp[— (2/a)| ln(r/a)|] and the corresponding vacuum stresses are strongly 
anisotropic: Ph ~ ap±h. Note that this suppression effect also takes place in the scalar 
case [9, 10]. 
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A Summation Formula over the Zeros of a Combina- 
tion of the Bessel Functions 



We have seen that the vacuum expectation values for the energy-momentum tensor for 
a spinor field inside a spherical shell on background of the global monopole spacetime 
contain sums over the zeros of the function 



X{z) = zJliz) +(^11 + S^^Z^+H^ - (-l)^i/) J^{z) , (7 = 0, 1 , 



±1 



(72) 



To obtain a summation formula over these zeros, we use here the generalized Abel-Plana 
formula [11, 12]. In this formula, as a function g{z) let us choose 



(73) 



with a function f{z) analytic in the right half-plane Re^; > 0, and Yij{z) is the Neumann 
function. For the sum and difference of the functions f{z) and g{z) one obtains 



/(;.) - {-lfg{z) 



f{z) , ^ = 1, 2 



(74) 



with h'^^ and H^^ being Bessel functions of the third kind or Hankel functions. By using 
the asymptotic formulae for the Bessel functions for large values of the argument, the 
conditions for the generalized Abel-Plana formula can be written in terms of the function 
f[z) as follows: 

\f{z)\<e{x)e'^y^ , z = x + iy, |^| ^ oo , (75) 

where c < 2 and e{x) for a; — >• oo. 

Let Xi,^s 7^ 0, s = 1, 2, 3 ... be zeros for the function Ju{z) in the right half-plane, 
arranged in ascending order, X^^^s < Xi,^s+i- By using the Wronskian 1^[J,^(2;), 1^1,(2;)] = 

2/7rz, one can easily see that Y{Xi,^s) = 2/ (^7rJ{Xu,s)J ■ This allows to present the residue 

term coming from the poles of the function g(z) in the form 



niRes^=x^^^9{z) = T^{Xu,s)f{Xu,s) 
where we have introduced the notation 

n{z) = — 



(76) 
(77) 



Substituting (73) and (74) into the generalized Abel-Plana formula [11, 12] and taking in 
this formula the limit a ^ (the branch points z — are avoided by semicircles of 
small radius), we obtain that for the function f{z) analytic in the half-plane Re^; > and 
satisfying condition (75) the following formula takes place 



{n „b 
y]T.(A.,,)/(A.,,)- / f{x)d: 

-"""'fixe 



I&\x) 

li-'Hx) 



+ e'^'^7(xe 



-7ri/2' 



Ki-\ 



X) 



&\x) 



dx 



(78) 
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where on the left X^^n < b < Xu,n+i, and Ti,{X^^k) is determined by relation (77). In formula 
(78) we use the notation 

^ " j zF'{z) + (/X ± sj^^^^ - i-iyiy) F{z) , \z\>ii 
for a given function F{z), and have taken into account that for z — iy 

( (/i^-y^)^/^ \y\<i^ 

(^2 + ^2y/2 _ J ^y2_ ^2y/2^^n/2 (gQ) 

[ (|/2-;x2)V2e-W2, y^_^ 

In this paper, we apply formula (78) to the sums over s in the expressions (35) for the 
vacuum expectation values of the energy density and vacuum stresses. As it can be seen 
from expressions (36)-(38), the corresponding functions f{z) satisfy the relation 

e-^'^7(xe'^^/2) = -e^'^7(xe-'^*/2) , for 0<x<n . (81) 

By taking into account that for these x one has = we conclude that the part 
of the integral on the right of Eq. (78) over the interval (0, //) vanishes. 
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